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ABSTRACT 


Applied Mathematics have been successfully used in the development of 
science and technology in 20th -21st century. In Mechanical Engineers, an 
application of Mathematics gives mechanical engineers convenient access 
to the essential problem solving tools that they use. In this paper, we will 
discuss some examples of applications of mathematics in Mechanical 
Engineering. We conclude that the role of mathematics in engineering 
remains a vital problem, and find out that mathematics should be a 
fundamental concern in the design and practice of engineering. 
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INTRODUCTION 

In this paper, several examples of applications of 
mathematics in mechanical engineering are discussed. 
Mathematics occupies a unique role in the Mechanical 
Engineering and represents a strategic key in the 
development of the technology. In this paper we elaborate 
some topics such as Matrices, Laplace transform, Partial 
differential equation for Mechanical Engineering. 


SOME OF THE MATHEMATICAL TOOLS THAT ARE 
USED IN MECHANICAL ENGINEERING 

> Matrices, 

> Laplace transform 

> Partial differential equation 


Matrices 

Matrices: A rectangular arrangement of number, symbols, 
or expressions in rows and columns is known as matrix. 
Matrices play an important role in mechanical engineering 
syllabus. Some subjects mention below in which we will 
apply matrix knowledge - 


In Engineering Material Sciences (Miller indices) matrix 
play an important role for defining crystal lattice 
geometries. In Strength of materials Strain matrix, stress 
matrix and the moment of inertia tensors are used for 
solving problems. We will also find application of matrices 
in analogous subjects like Design of Machine Elements, 
Design of Mechanical Systems. MATLAB stands for “Matrix 
Laboratory’. Matrix is basic building block of MATLAB. 
Computer -aided Designing (CAD) cannot exist without 
matrices. In Robotics Engineering it is impossible to 
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design a robot without the use of matrices. All the joint 
variables for forward/inverse kinematics and dynamics 
problems of the subject are noting down by matrix. 
Similarly, many concepts of matrices are used in Finite 
Element Analysis (FEA) and Finite Element Methods 
(FEM) for solving problems, just like CAD does. Mainly 
eigen value concept of matrices is used here. 


These are some applications of matrices in mechanical 
engineering. Now we discuss a one example of matrix 
representing stress for calculation of principal stress. 


INTRODUCTION TO THE STRESS TENSOR 
Y 
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Tyx Oyy Tyz 


s Txy rn 
O7z 


Ux Tey 
Here diagonal elements are normal stresses and off 
diagonal elements are shear stresses. On stress element 
in 3D, the normal and shear stresses can be compiled into 
a3x3 matrix is  calledastress tensor. From_ our 
observation, it is possible to find a set of three 
principal stresse for a given system. We know that, the 
Shear stresses always become zero when principal 
stresses are acting. So in view of stress tensor, according 
to mathematical terms, this is the © process 
of digonalization of matrix in which the eigenvalue of 
given matrix play the role of principal stresses. 


Example 

The state of plane stress at a point is represented by the 
stress element below. Find the principal stresses and 
angles at which the principal stresses act. 


—_ 





uetA=[n° og l=Lt0 a2! 


Cy 12 
Now we calculate the eigen value of matrix A. 


We consider the matrix form AX=AX, where I is identity 
matrix and A become eigen value, X is eigenvector such 
that (A-AI)X=0. 


To find the value of A we consider Characteristic equation 
which is given by det(A-AI)=0 
(—20 — A) —10 | _ 0 
—10 (12-—A)| © 
A’ + 8-340 =0 
A = 14.8679, —22.8679 


So the principal stresses are 14.8679 and -22.8679 as 
discuss above. By using eigen values, we can calculate 
eigenvectors. These eigen vector are use for calculating 
the angle at which the principal stresses act. 


Let X= | be the eigen vector such that (A-AI)X=0 
i — A) —10 | ] _ °) 

—10 (12 —A)} l%2! Lo 
For 


A = 14.8679 oe 


0 


— 14,8679 ) ~10 | ia _ 
40 (12 — 14.8679 )| lx2] ~ 


a rem | Fe ~ H 


On solving above matrix, first eigenvector is given by 
[9.2807] 
1 


For 
A= 
7 (—20 + 22.8679 ) —10 x4] 
eo —10 (12 + 22.8679 ) | 7 
i 
0 
(2.8679 ) —10 


—10 (34.8679 1 Fe ~ i 


On solving above matrix, second eigenvector is given by 
nea 
1 


Now we calculate angles at which the principal stresses 
act, but before that we can check whether eigenvectors are 
correct or not. For this we will digonalized given matrix as 
fallow 


” (0.265 oa |< a | oo taal 
0.265 0.07591'-10 12 1 1 


; oun L=4 


Here B is model matrix (matrix combination of all 
eigenvector) and D is digonalized form of matrix A 


These imply that the calculated eigenvector are correct. 


Now to calculate angle, we must calculate unit eigenvector. 
wowed S$ vohed and aa e eae 
1 0.9612 1 0.2756 


Now compile this unit eigenvector into rotation matrix R 
such that determinant R = +1 
R= 0.9612 ae | 

0.2756 0.9612 


Determinant R = (0.9612) (0.9612) + (0.2756) (0.2756) = 
1 


_[cos@ —sin@ 

R= ond eal 

D’=RTAR 

_ 0.9612 oe 7 | 0.9612 Pia 
—0.2756 0.9612!t—-10 12 110.2756 0.9612 


~ ou ren 


So 8 = 16°, as we found earlier for one of the principal 
angles. Using the rotation angle of 16°, the matrix A 
(representing the original stress state of the element) can 
be transformed to matrix D’ (representing the principal 
stress state). 
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20MPz 





\ 


In this way we can calculate the principal stresses and 
angles at which the principal stresses act by matrix 
method. 


Partial differential equations 

A large number of problems in fluid mechanics, solid 

mechanics, heat transfer, electromagnetic theory and 

other areas of physics and engineering science are 

modelled as Initial Value Problems and boundary value 

problems consisting of partial differential equations. In 

this paper, some of most important partial differential 

equations of one dimensional heat equation have been 

derived and solved. 

> Partial differential equations are used for heat 
conduction analysis. 

> Second order differential equation is used to find 
maxima and minima of function of several variables. 

> Partial differential equation help to provide shape and 
interior, exterior design of machine. 

> Partial differential equation is used to calculate heat 
flow in one and two dimensions. 


One dimensional heat flow 

Let us consider a conduction of heat along a bar whose 
both sides are insulated. Also the loss of heat from the 
sides of the bar by conduction or radiation is negligible. 
One end of the bar is taken as origin and direction of heat 
flow is along positive x-axis. The temperature u at any 
point of the bar is depend upon the distance x of the point 
from one end and the time t. The temperature of all points 
of any cross-section is the same. 





Hence, the quantity of heat Q: flowing into the section at a 
distance x will be 


Ou 
QO; = -Ka{ 


The negative sign on RHS is because u decreases as x 
increases, 


x Per second 


The quantity of heat Q2 flowing out of the section at a 
distance * + 9X will be 


Ou 
Q, _ -1a{ 


x+ox Per second 


The amount of heat retained by the slab with thickness ox 
is, therefore, 


0,-0, = KA (Se) ; er), 


Per second 
Geet ete aiicttee oocenees (1) 
The rate of increase of heat in the slab is 
Ou 
= sopAox—, 
OF cca psetscieeauetn anions (2) 


Where s is the specific heat and p is the density of the 
material of the bar. 


From (1) and (2) we have therefore 
7) a a 
spades = KAL(s.) soy ~ (onde 


[aac 


ot So 6x 

du 0*u K 
= 6.072. where C4, = — 
Ot ax2 ’ Sp 


is known as the thermal diffusivity of the material of the 
bar. 


Solution of heat Equation 
The heat equation is 


du , 0*u 
Ot 0x? 
where the symbols have got their usual meanings. 
Let uw = T(t)X (x), oe cee cee eee 2 
Then = x2 
at dt 
And 
d*u  d*x 
ax? dx? 


Taking the above substitutions in (1), we obtain 
dT d?X 


X= CT 
x 


1 d*X K(say) : 
Sse = Se SY) srr: 

CT dt Xdx2 9° > 
The solutions will now be found under the following three 
cases: 
Case I: When K=0, we have from (3) 


1 4 1d*X _ | 
C?Tdt = 'Xdx* — 
aT 9 UX , 
> ——_ = 
dt "dx? 





@IJTSRD | Unique Paper ID - JTSRD38348 _ | 


Volume-5|Issue-2 | 


January-February 2021 Page 144 


International Journal of Trend in Scientific Research and Development (IJTSRD) @ www.ijtsrd.com eISSN: 2456-6470 


After integration we get 
= T = C,and X = (2x + C3 


Using this in equation (2), we get 
Which is a solution of (1). 


Case II: When K= m4, i.e. K is >0, we have from equation 


(3) 

1 dT  ,1d*x_ , 
C2T dt XG? = 7 

d d*X 

ee 272 ich 2 = 
=> r m‘*C~dt, dx? m 0 


= logT = m*C*t + logC, ,A.E.D* — m? = 0 


= T= Ce™t X =C.F.= Coe™ + Cge7™ (~ Pl. 
— 0) 


From (2), we have therefore 
>u= Cye™Ctese™™ tego) OB) 


which is a solution of equation (1) 
Case III: When K= - m2, i.e. K < 0, we have from (3) 
1 dT 5.1 d?X . 


Crd" Xde ~~ 
aT aX 
op = Tm Cat, a tmx = 0 
=> logT = —m’*C?t + logC,, 


AE. D*?+m* =0>D=+4mi 
=> T = CreTm ct 


=“ fi => a= 0,f =m 


X= C.F.=e** (Cg cosBx + CosinBx) 
=e°* (Cg cosmx + Cysinmx) 
= Cg cosmx + Cosinmx (° P.I.= 0) 


From (2) we have 
u= Crem tC, COS MX + Co SIN MX) re see ee ee ees (C) 
which is a solution of equation (1). 


Among these solutions, we have to choose that solution 
which is consistent with physical nature of problem. Since 
u decreases as t increases, the only suitable solution of 
heat equation (1) is solution (C). 


Example 

Example: Find the temperature in a bar of length 2 units 
whose ends are kept at zero temperature and lateral 
surface insulated if the initial temperature is sin + 


3 sin = 
Solution: One dimensional heat equation is 
du 02 07u , 

ap © aye (1) 


The solution of equation (1) consistent with physical 
nature of problem is given by 


OCT) SS0Gl 0 case ees) 
UWOGD) = 0G XS 2 oes itcces 4) 
1X 51x 
u(x,t) = sin > + 3 sin —— Ott = Oesuieanlo) 


Using the condition (3) in (2), we obtain 
0= Cem (C,) 
=> C, = 


From (2), we have therefore 
U = CyeT™ EC, SIM MX eae vee ses vee vee one (6) 


Using condition (4) in (2), we obtain 
0=Ce"™ "tC, sin 2m 

=> sin2m = 0 = sinnt 

=> 2m = nt 

=>m=nt/2 


From (6) we have 
nly? 224 NIX 
_ -(+-) c 
u=Cj\e (=) GC sin(— , ) ee se 


The general solution is 
00 
~ - (28) ct _ (ntx 
u= » bye \2 sin (—) ieee eee CO) 
n=1 


Using (5) in (8), we obtain 
oe) 


TUX _ 5mx (NIX 
sin +3 sine = yy b, sin (—) 


n=1 
is b, sin (= =) + by sin (= =) +b; sin (= a) + by sin (=) + 
ment 1 =p a) = b, =b, = be =~ = 0 


From (8), we have therefore 


Rigs, Laipx ~ (82) ce 57x 
u=b,e 4° sin—-+bDese \2 sin —— 
7 2 
2 2 
19 TUX STL)" 24 510X 
Meo “t . =(=-) C 
Ty e>-4 sin + 3e (=) * sin—— 


which is the required temperature. 


Laplace Transform 

Laplace Transform: The Laplace Transform is_ the 
transform to time domain (t) into complex domain (s). 
Laplace Transform plays an important role in engineering 
system. The concept of Laplace Transform is applied in the 
area of science and technology such as to find the transfer 
function in mechanical system. 


The Laplace Transform of f(t) is defined and denoted as 
[FO] =f, eo fO at 
Where f(t) is the function of t and t > 0. Provided that 


integral exist, s is the parameter which may be real and 
complex. 





uUu= Crest (C, cos mx + C3 SiN MX) ... ee (2) Find the transfer function of the system shown below? 
Where 
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Transfer function: It is the relation between the output 
and the input of a dynamic system written in complex 
form (s). For a dynamic system with an input u(t) and an 
output y(t) ,the transfer function H(s) is the ratio between 
the complex representation (s) of the output Y(s) and 
input U(s). 


Whenever they give any mechanical translation system, 
mass desk pot, spring will be in mechanical system. When 
we have to find out the transfer function of the system we 
need to take the output transfer by input transfer .Output 
abst 

1(s) 


is in terms of x, and input in terms f,. If we consider 
we get transfer function of the system as fallow 


First mass m,, second mass m, 


a? x4 
dt? 





Force due to the mass m, is fm, = m, 


Taking Laplace Transform both side 


Force due to the spring k, is fk, = k,x, 
Force due to the desk pot D, is fb, = b, < (xX, — X2) 
Force due to the spring k, is fky = ky(x, — X2) 


External force is equal to the sum of all internal forces. 
fit) = fm, ue + fb, + fkz 

d d 
fii@© = mM > + kx, + bs (X1 — X2) + kg (x1 — x2) 


Taking Laplace transform both side, we get 

Fi(s)= m,s* X,(s)+ ky X, (s) + by s CX, (5s) - 
X, (s) ) + kz CX; GS) — X2 (S)) 

Fi(s) =( m,s* + biyst ky + k2)X,(s) - (bbs + 
kz) X> (s). 

Force acting on mass m,. 


d*x> 
dt? 





Force due to the mass m, is fm, = my, 


Force due to the desk pot D, is fb, = b, < (X2 — X) 


ax2 


Force due to the desk pot bj is f bz i 


Force due to the spring k, is fk, = k2(x2 — x1) 
No external force acting on mass m, 
0=fm,+fb,+ fb, + fk, 


d*x d dx 
as + b, Ae (x2 — x1) + by + Ko (x2 — x4) 





0 =m, 


0 = Mz S* Xz (s) + b, s (X; (s) — Xy (s)) + by s Xz (Ss) + kz (X2 (S) — X; (S)) 


0 — X> (s) (m,s7 +b, st+b,s + ky ) = X,(s) (bi s+ k, ) 


(mz s* +byS+b25+kz ) X2(s) 
(by s+ kz) 


X, (s) 


Fi (s) =(m,s*+b,st+k,+k,) 





(mz s* +byS+b25+k2) X2(s) 


— (bys + kz) Xz (Ss) 


(by s+ kz) 
_ (m, 5 + by S+ky+ K2 ) (m2 s? +b,S+b2s5+k2 )X2(s)- (by St+ kp) 
Fy (s) = 4218) 9) OOOO Oe 
(by $+ kz) 
X2(s) _ { (b1 $+ Kz) 
Fy (s) (my, s2 + by St+ Ky+ K2 ) (m2 s2 + by S+t+ b2 S + K2 )X2(s)- (by S+ kz)? 


It is a transfer function of given mechanical system. 


CONCLUSION 

In this paper we conclude that mathematics is backbone in 
study of technical subject of Mechanical Engineering. 
Mathematics is applied in various field of mechanical 
engineering like maths is use in fluid mechanics, straight 
of material, machine design etc. In this way mathematical 
concept and procedure are used to solve problem in above 
mentioned field. 
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